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Abstract 



In the Sobolev space ^ (0,1) optimal quadrature formulas with the nodes (1.5) are inves- 
tigated. For optimal coefficients explicit form are obtained and norm of the error functional is 
calculated. In particular, by choosing parameter 770 in (1.5) the optimal quadrature formulas with 
positive coefficients are obtained and compared with well known optimal formulas. 
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1 Introduction. Statement of the problem 

It is known, that numerical integration formulae, or quadrature formulae, are methods for the 
approximate evaluation of definite integrals. They are needed for the computation of those 
integrals for which either the antiderivative of the integrand cannot be expressed in terms of 
elementary functions or for which the integrand is available only at discrete points, for example 
from experimental data. In addition and even more important, quadrature formulae provide a 
basic and important tool for the numerical solution of differential and integral equations. 

There are various methods in the theory of quadrature, which allow us approximately cal- 
culate integrals with the help of finite number of values of integrand. The present paper also 
is devoted to one of such methods, i.e. to construction of optimal quadrature formulas for 
approximate evaluation of definite integrals in the space L% (0, 1) equipped with the norm 
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Consider a quadrature formula 



r N 

J tp{x)dx^Y. C W{xp) (1-1) 



(3=0 

with the error functional 

N 

£{x) = e m (x) - C P S ( X - x p) (1-2) 

(3=0 

in the Sobolev space L^ifi, 1). Here and a;^ are the coefficients and the nodes of the 
quadrature formula (1.1), respectively, S(x) is the Dirac's delta-function, £[ 0) i](x) is the indicator 
of the interval [0,1]. 

In order that the error functional (1.2) is defined on the space L^(0, 1) it is necessary 
performing following conditions (see [28]) 

(£(x),x a ) = 0, a = 0,1,2,..., to- 1. (1.3) 

The difference 

1 N oo 

(£(x),<p(x))= ip(x)dx-J2 c M x p)= / t(x)<p(x)dx. (1.4) 
o 

is called by £/ie error of the quadrature formula (1.1). 
By Cauchy-Schwartz inequality 

| (€(*), V (x))| < |b(x)|^ m) (0, 1)|| • ||€(x)|4 m) *(0, 1)|| 
the error (1.4) of the formula (1.1) is estimated with the help of norm 

*(s)l4 m) * = sup |(€(x),^(x))|. 

1^)14™) ||=i 

of the error functional (1.2). Consequently, estimation of the error (1.4) of the quadrature 
formula (1.1) on functions of the space ^^(0,1) is reduced to finding of norm of the error 
functional £(x) in the conjugate space L^^O, 1). 

Clearly, that norm of the error functional £(x) depends on the coefficients Cp and the nodes 
Xp. The problem of finding of minimum of norm of the error functional £{x) by coefficients Cp 



and by nodes xp, is called by S.M.Nikolskii problem, and obtained formula is called optimal 
quadrature formula in the sense of Nikolskii. This problem first considered by S.M.Nikolskii 
[15], and continued by many authors, see e.g. [2-4,16,32] and references therein. Minimization 
of norm of the error functional £(x) by coefficients Cp when the nodes are fixed is called Sard's 
problem. And obtained formula is called optimal quadrature formula in the sense of Sard. First 
this problem investigated by A.Sard [17]. 

There are several methods of construction of optimal quadrature formulas in the sense of 
Sard such as spline method, (p— function method (see e.g. [23], [2]) and Sobolev's method which 
based on construction of discrete analogue of a linear differential operator (see e.g. [26]). In 
the different spaces, based on these methods, the Sard's problem investigated by many authors, 
see, for example, [1,2,4,6,7,10-14,19-26,28-31] and references therein. 

Furthermore, in works [11,19,21,31] were obtained explicit formulas for coefficients of optimal 
quadrature formulas for any m and for any number N of the nodes Xp in the space 

By I. J.Schoenberg and S.D.Silliman [25] investigated the Sard's problem for the case N — > oo 
in the space L^. In [25] an algorithm for finding of optimal coefficients was given with the 
help of spline of degree 2m — 1. In the cases m = 2,3, ...,7 the coefficients are calculated 
using a Computer. There was especially noted, that when m = 7 among optimal coefficients 

(7) 

the negative coefficient B\ appears. It appears that by increasing m the number of negative 
coefficients increases. This is confirmed by computations of the optimal coefficients for m < 30 
in [30]. 

But, as known that in applications the optimal quadrature formulas with positive coefficients 
play very important role. 

Naturally, it is arises a question: can we obtain with some way the optimal quadrature 
formulas with positive coefficients in the Sobolev space 1,^(0, 1)? 

The main objective of the present paper is, using the Sobolev's method, to construct of 
optimal quadrature formulas in the sense of Sard in the space L^^O, 1) with the nodes 

Xi = r)ih, x N _i = 1 - rjih, i = 0, t - 1, < r] < rji < ... < r) t -i, 
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xp = h(3, t<(3<N-t, h 



1 

N' 



Y when m is even, 

[y] + 1 when m is odd, 



(1.5) 



where [a] is integer part of the number a and choosing the parameters rji, i — 0, t — 1 to obtain 
the optimal quadrature formulas of the form (1.1) with positive coefficients. This means to find 
the coefficients Cp which satisfy following equality 



e(x)\L 



(m)* 
2 



= inf 



e(x)\L 



(m)* 
2 



;i.6) 



for the nodes (1.5). 

Thus, in order to construct of optimal quadrature formula in the sense of Sard with nodes 
(1.5) in the space L^iO, 1) we need consequently to solve following problems. 

Problem 1. Find norm of the error functional £(x) of quadrature formulas of the form 
(1.1) in the space L^*(0,1). 

Problem 2. Find the coefficients Cp which satisfy the equality (1.6) with nodes (1.5). 

It is known [26,28,29], that square of the norm of the error functional (1.2) for arbitrary 
fixed Xp has following form 

2 



x)\L^*(0,l) 



N N 



K f3=0 7=0 





|2m— 1 


\ x - 





2- (2m- 1)! 



N 
13=0 



X 



2m 



+ (l-x f 



\2m 



+ 



(1.7) 



2 -(2m)! (2m + 1)! 

Obviously, that the square of the norm (1.7) of the error functional £(x) is multidimensional 
function with respect to the coefficients Cp . Moreover, the error functional £(x) satisfies the 
conditions (1.3). In [28,29] taking into account these facts the Lagrange function with conditions 
(1.3) is constructed for finding condition minimum of the (1.7) and differentiating that function 
by Cp and by \ a (where \ a are Lagrange factors) following system of linear equations was 
obtained 



N 



E \Xp Xj 
7 2 • (2m 



2m-l 



m—1 



7=0 



1 



X 



2m 



+ (l-xp) 



2 m 



o=0 



0! 

N 
7=0 



2 (2m)! 



xp e [0,1], 



a 



a = 0, m — 1. 



(1.9) 



It was proved in [28,29] that this system has unique solution and this solution gives minimum 
to the expression (1.7). This means that square of norm of the error functional £(x) being 

o 

quadratic function of the coefficients Cp has unique minimum in concrete value of Cp =Cp- 

o 

The quadrature formula with the coefficients Cp when the nodes are fixed is called optimal 

o 

quadrature formula in the sense of Sard and the coefficients Cp are called optimal one. 

o 

Below for convenience the optimal coefficients Cp we shall remain as Cp. 

Thus, the problem 1 is already solved by S.L.Sobolev in the space and the problem 2 
is reduced to the system of linear equations for optimal coefficients. Should be noted, that in 
[28,29] the problem 1 solved for multidimensional case, i.e. for cubature formulas. 

In the present paper we will solve the system (1.8)-(1.9), i.e. we will solve the problem 2. 



2 Definitions and known formulas 

In this section we give some definitions and formulas which are necessary in the proof of the 
main results. 

The Euler-Frobenius polynomials E k (x) , k — 1,2, ... are defined by following formula (see, 
e.g. [26]) 



1 - x) k+2 ( d 



k 



x 



E„{x) = 1. 

For Euler-Frobenius polynomials following identity holds 



E k (x) = x k E k ( - ) , (2.2) 



and also following theorem is take placed 

Theorem 2.1 [20]. Polynomial P k [x) which is determined by formula 

ft(x) = («-l)W^^ 15f (2.3) 

is the Euler-Frobenius polynomial (2.1) of degree k, i.e. P k (x) = E k {x). 



5 



Following formula is valid [9]: 

n - 1 I k / _ v « n fc / « X* 



I7=n, 



(2.4) 



7=0 J i=0 v J/ a i=0 

where A*7 fc is finite difference of order i of 7 fe , A l fc = A*7 fc | 7=0 . 

At last we give following well known formulas from [8] 

p-i fc+i 



7=0 j=l J v Jy 



where B^+i-j are Bernoulli numbers, 



A a ^ = ^ j ^ I A a p x u - p . (2.6) 

3 Auxiliary results 
3.1 Lemmas 

In the proofs of the main results we need following lemmas. 

Lemma 3.1. [20] The Optimal coefficients of quadrature formulas of the form (1.1) with 
the error functional (1.2) and the nodes (1.5) in space L 2 "^(0, 1) have following form 

(rn-l \ 
i+E^+^n (3.i) 

when (3 — t,t + l, N — t. Here dk are unknown parameters, qu are roots of the Euler-Frobenius 
polynomial E2 m -2{q) , which \qk\ < 1. 
Lemma 3.2. Following identity holds 

y tiiLgfc A*0" = (-1)" V Qk +l lj * k A*0 a , (3.2) 

where a = 1,2, m — 1, s = 0, 1, gfc are roots o/ £/ie Euler-Frobenius polynomial £^-2(0)- 
Lemma 3.3. T/ie coefficients of optimal quadrature formulas of the form (1.1) with the 
error functional (1.2) and the nodes (1.5) in the space L^^O, 1) satisfy following system 

*-! Z*- 1 na a ( _-i \i„t+i _ „N-t+l \ 

E cw? = a E /?" + t + E * E ' («> 

/9=0 \/3=l fe=l i=0 Kqk ' J 



, 0, 2, 4, m — 2 when m — even, J 1, a — 0, 

/iere ct = < CT = < 

0, 2, 4, m — 1 when m — odd, I 0, a^O, 

, Y w/ien m even, 

£ = ^ where [a] is integer part of the number a, c4 are unknown 

[y] + 1 w/ien m is odd, 

parameters, are roots of the Euler-Frobenius polynomial E 2m -2(q) , \Qk\ < 1- 

Lemma 3.4. Let m > n and m,n & N then following formula for binomial coefficients is 
true 

Cm—l/~m sim—2/~m , , / -i \m— n— 2/~m+l/~m , / i\m-ti-l/TO/in /~m (1 A\ 

m Vl ~~ U m L/ m-2 + - + r i J °m U n+l + l _i J U m°n ~ °m' V 6 -^) 

Lemma 3.5. [28] Following identity is true 

y k = lt 7r Al ° fc > (3 - 5) 
i=i 

w /iere yW = y(y - l)(y - 2)...(y -k + 1), A*0 fe = £ (-l^-'CjZ*. 

i=i 

3.2 Proofs of lemmas 

Lemma 3.1 is proved in [20]. 

Proof of lemma 3.2. We will denote the left and the right hand sides of the identity (3.2) 
by A\ and A 2 , respectively, i.e. 

01 n N ~ s+1 4- ( 1 Yn s+i a a s+1 4- ( 1 Y+ 1 n N ~ s+i 

Al = y z3h +l-_V gfc A 'o Q , a 2 = y Qk +{ ij .^ fc A*o a . 

Then, using theorem 2.1 and the identity (2.2), for A 1 we obtain 



_ q N - S+ l + ( _ 1)agj H-l 



= ** ; ^ ga-i(ft) • (3.6) 



For A 2 also, using theorem 2.1 and the identity (2.2) we have 

_ q N -s +i + ( _ 1)ttgj H-l 



^ = Vfc ; Aa+l ^ ^a-l((ft) • (3.7) 



From equalities (3.6) and (3.7) we will obtain the statement of lemma 3.2. Lemma 3.2 is proved. 



Proof of lemma 3.3. Here we will use lemma 3.1 and the equality (1.9). The optimal 
coefficients Cp and the nodes xp have to satisfy the equality (1.9), i.e. 

N 1 

E^^q— I' « = 0,l,2,...,m-l. (3.8) 

From symmetry of the nodes (1.5) follows that 

Cp = C N .p. (3.9) 

Here the case a = we consider separately. 
Let a = then from (3.8) we have 

N 

$^ = 1. (3.10) 

p=o 

Taking into account (3.9) and using lemma 3.1 for the left side of (3.10) we get 

N t-1 N-t / m-1 \ 

P=0 p=o p=t \ fc=l / 

t- 1 ( m-1 t _ n N-t+l\ 

Hence, taking into account (3.10), we obtain 

E^^f^-Ei^fV ' (3 ' n) 

/3=0 V fc=l yfc / 

This is the first equation of the system (3.3) corresponding to the case a — 0. Thus we have 
proved the lemma for the case a — 0. 

Let, now, a = 1,2, m — 1. Then, using symmetry of the nodes (1.5) and keeping in mind 
the equality (3.9), for the left side of (3.8) we obtain 

TV t-1 N-t 

E c ^ = E °p M + ( x - + E c ^ = Y ± + y * 

/3=0 /3=0 /3=t 

where 

t-l N-t 

Y x = EQ + (1 - ^)°) , E 2 = EW 

/9=0 /3=t 



Substituting the expression (3.1) of the optimal coefficients Cp into Y 2 , adding and subtracting 
the expressions 

t-i / m-l \ 

EM 1 + E^(rf+^) (W 



and 



we have 



/3=1 
AT-1 



fc=l 



m— 1 



E Mi+E*(s!+^) <w. 

=AT-t+l V fe=l / 



p=N-t+l 
<N-1 



m-l AT-1 

y 2 = h~» ( E 0° + E M E ^ a + E ^ 

\/3=l fe=l \/3=l /3=1 

t-1 / m-l 

-aE «w + (i - i + E d * (rf + ^) 



0=1 



fc=l 



Using equalities (2.4), (2.5), (2.6) and binomial formula, then grouping in powers of h and 
taking into account lemma 3.2, for Y% we have 



Y 2 = h a+1 



m— 1 a 
k=l i=l 

Q-l 



(l"9 fc ) 



t-i 



m— 1 



i+1 



k=l 



a^-e^(i + E^ (rf+^-'O 

1)^ 



/3=1 
m-l i 



(i+(-in+ 



i + E^E 



+ (-iy +l E/*(i + E d *(rf+flO 



t-i 



0=1 



fc=l 



2t- 1 ^,4-?r m ' 



-E^ 



fe=i 



+ 



a + 1 



Now we consider Y\. Taking into account that when (3 = 0,t — 1 then from (1.5) the nodes 
X/3 = rjph and applying binomial formula, after that, grouping in powers of h, for Y\ we have 



t-i 



t-i 



y, = E W + (! - = E M + E 



0=0 



0=0 



t-1 



a-1 



= /i-+ i Efc- 1 c^(i+(-i) a )+E' 



W +1 



3=0 

t-1 



j\(a — j)\ 

0=o j=o J v J ' 0=0 
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Hence using (3.11) we obtain 



t-i 



0=0 



+ E 



2«-l ^-ri-tf-* 1 ' 



E^ 



fe=i 



Now, adding Fx, F 2 and the result substituting to the left side of the equality (3.8), we obtain 

m— 1 a 



r+ i 



E*E 

fc=i i=i 



(1 - <7*) i+1 



AW 



t-l / m-l \ t-1 

E ? a 1 + E d * (rf + ^"") + E * _1 <w 



0=1 



fc=l 



t-i 



5 



/3=0 

m-l j M +i 



(l + (-l) Q ) + 



J + 
m— 1 



^ + E*E 



fc=i i=i 



+ (-l)Vft 



A l J - 



t-1 / m-l \ t-l 

E W 1 + E * (*£ + + E ^(-^) J 

0=1 V fc=l / 0=0 



0. 



(3.12) 



Clearly, that the left side of (3.12) is the polynomial of degree a + 1 with respect to h. Since 
this polynomial is equal to zero then all coefficients of the polynomial are zero. 
When a is odd from (3.12) we get system of equations 



m— 1 



t-l 



fc=i 



i=l 



5 



t-l 



'J'+l 



i + (-1)' E (C/j/i- 1 ^ - P j ) =0, j = 1, 2, a - 1, 

/3=0 



which is the part of the system (4.22), because 1 < a < m — 1. 
When a is even from (3.12) we get following system of equations 

m_1 f j n N + i A_ ( iVn * _1 / N 1 

E * E = 7i + s> (rf + «r") + 

fc=l I i=l ^ W 0=1 J 



5 



t-i 



0=0 
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which is the part of the system (4.22) and one new equation for each even a 



m—l 

E* 

fc=l 



a N+i | / -. \ j t—1 

i=l 1 0=1 



t-1 



0=0 



Hence and from (3.11) for the optimal coefficients Cp, (3 = 0, ...,t — 1 we will get following 
system of equations 



t-i 



m— 1 



E^l^-E*^ 



t-i 

/3=0 



0=0 
m—l 



k=l 



„N-t+l 
Qk 

t-1 



t—1 m—l / a N+j I / -1 \i * — 1 

5> - £ * E ~\ - E ^ (rf + 

0=1 fc=l \i=l V yfe; 0=1 



N-13 



where a 



0=1 fe=i 

2,4, ...,m — 2 when m — even, 
2,4, ...,m— 1 when m — odd. 



(3.13) 
(3.14) 



Thus for solvability of the system (3.13), (3.14) (with respect to (Tg, /3 = 0,t— 1) t have to 
take on a value y when m is even and [y] + 1 when m is odd, where [a] is the integer part of 
the number a. 



t-i 



Now, let Y 3 = J2P a [Qk + Qk'^j- Th en, using formula (2.4) and lemma 3.2 taking into 
/3=i ^ ' 
account that a is even natural number, we have 



-Qk + 9* A in« , y~ L > % -Ik 



3 4^ fl-g fc ) i+1 Z - 



(3.15) 



i=l v Jn,/ i=0 

Putting the equality (3.15) into (3.14), taking into account (3.13) we will get the assertion of 
lemma 3.3. Lemma 3.3 is proved. 



Proof of lemma 4. Using formula of the binomial coefficients we will get 

m 1 ~ l\-(m-l)\k\-(l-k)\ ~ 
~ ml im ' k)l ~C m C^ k . (3.16) 



k\-(m-k)\ (I - k)\ ■ ((m — k) — (I — k))\ 
Applying the equality (3.16) to the left side of (3.4) we obtain 



^jm—l^jn ^jm—2^jn _|_ _|_ / j\m-Ji-2^m+l^in _|_ / |\m-n-l^n^n 
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^m^rrt— ra ^m^m— n ~r ••• ~r \ -LJ ^m^m— ra 

Cn I r~im— n— 1 sim—n—2 , , / i \m- rt— 1/^rO \ 

m l°m-n ~ ^m-n "+" ••■ "+" l — -U ^m-nj ~ 

Cra I rtra—n ( rim~n r~im—n—l i i / i \m-ii/n0 \\ 

m \ym-n ~ \ym-n ~ ^m-n "+" "■ "+" l — X J U m-njJ ~ 

Cn I r~im—n f-i i \m-n\ /-in 
m \ym-n ~ K 1 ~ l ) ) ~ U m) 

Lemma 3.4 is proved. 

Proof of lemma 3.5. The assertion of this lemma we will prove by induction. 



For k — 1: 



v [i\ 7/ [l] 

»=El A ' ' =^'0'=!,, 



1=1 

For fc = 2: 

^ 7/M 7/W 7/ PI 

i=i 

Now, suppose that the equality (3.5) is true for k = a, i.e. 

a = y- y^ * = y A i « + *_A 2 a + ... + — — — A a a . (3.17) 
y ^ i\ l! 2! a! v 7 

i=l 

We will show that the equality (3.5) is also true for the case k — a + 1. For this we will use 
following equalities 

y^=yW(y-a), (3.18) 
A*0 Q+1 = i (A i_1 a + A^O") . (3.19) 
Multiplying both parts of the equality (3.17) by y and keeping in mind (3.18), (3.19), we obtain 

£±i 7/H 7,W 7, PI Tl [a+1] 

a+i _ £_A*0 a+1 =— A 1 Q+1 + — A 2 a+1 + + — A a+1 Q+1 

y ^ i\ 1! 2! "' (a + 1)! 

i=i v ' 

Thus, for any k E N the equality (3.5) is true. Lemma 3.5 is proved. 



12 



4 The main results 



4.1 Coefficients of optimal quadrature formulas 

For the coefficients of optimal quadrature formulas of the form (1.1) following theorem holds. 

Theorem 4.1. The coefficients of optimal quadrature formulas of the form (1.1) with the 
error functional (1.2) and the nodes (1.5) in Sobolev space £^(0, 1) are expressed by formula 

where dk, (k — l,m — 1) satisfy following system of m — 1 linear equations: 

m -! 3 _t+l , / 1 \i n N-t+i R <_1 

E*E ~* + ( ^ M = L_^±l _ V c k-i(t - mf , (4.2) 

fc=l i=l Wfc J 7 + 0=0 

j = l,2,3,...,m- 1, 

/jere the coefficients Cp = C^-p (/3 = 0, 1, ...,t — 1) are determined from the system (3.3), qu 
are roots of Euler-Frobenius polynomial E 2m -2(<l) , \qk\ < 1- 

Proof. As said before the optimal coefficients Cp {(3 = 0, N) are solution of the system 
(1.8)- (1.9). In the lemma 3.1 we have obtained representation of the optimal coefficients Cp for 
(3 = t, N — t. In the proof of the lemma 3.3, using the result of lemma 3.1, we have obtained the 
system (3.3) for the optimal coefficients Cp (@ — 0, t — 1). So we conclude that in order to solve 
the system (1.8)- (1.9) it is sufficient to find unknown parameters dk and unknown polynomial 

m— 1 

P m -i(%p) — X] ^aXp. For this in (1.8) instead of Cp, (3 = t, N — t substituting the expression 
fc=i 

(3.1) we will get polynomials of degree 2m with respect to xp in both sides of (1.8). Equating 
coefficients of same powers of xp we will get the system oi m — 1 linear equations for unknowns 
dk and we will find the coefficients of unknown polynomial P m -i(xp). Thus the proof will be 
complete. 

Further we will give detailed explanation of the proof of the theorem. 
Now we consider the equality (1.8) 

N |„ _„ \2m-l T 2m , (-i _ „ \2m 

E^^i)T + = 6 V ( 2 m y ™ <«) 
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We denote by 



\1m-\ 



yv P) ^7 2. (2m- 1) 

7=0 v 7 



/M = - 



2m 



2 (2m)! 

Now we consider the function g(xp) and let t < (3 < N — t. Then 



N 



gM = t c^r^T* E fi'T'C- E c 



\2m-l 



(2m- 1)! 

7=0 v ' 7=t 



(2m -1)! ^ 7 2(2m-l)! 



Further we denote by 



2m- 1 



(2m -1)! ' 



^2(^/3) = £c. 



7=0 

(^ - x,) 2 ™- 1 



7 =t 
N 



(2m -1)! ' 



*M = -E^%irr 



2m-l 



7=0 



2(2m- 1)! 



In the equality (4.8), using (3.1), adding and subtracting following expression 

2 m / 1— 1 m— 1 / 1— 1 t— 1 



/l 2 



(2m- 1)! 



£ (0 - 7) 2m " + £ M £ ^ - 7) 2 ™- 1 + £ - ^ 



.7=1 



we get 



h 2m 



(2m- 1)! 
/i 2m 



fe=i \ 7 =i 



7=1 



(4.4) 
(4.5) 

(4.6) 

(4.7) 
(4.8) 
(4.9) 



■ m-l / P P 

E (p - ^ + E 4* £ «i w - 7) 2 " 8 - 1 + E ^ w - 7) 2 " 1 " 1 

.7=1 fc=l \7=1 7=1 



t-1 m-l * _1 

E (/? - 7) 2 ™- 1 + E * £ «i - ^ + E (/? - ^) 2m_1 

.7=1 fc=l \7=1 7=1 



(2m - 1) 

Hence replacing (5 — 7 with 7 and using equalities (2.4), (2.5) 

h 2m 



^2{xp) 



(2m- 1)! 



V- (2m- l)!/j 2 m-j ^ 
^ j!(2m-j)! P 

_/3 2m-l 

-E 



m -l f / rl 2m "l Aj n 2m-1 1-/3 2m -l \io2m-l 
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2m- 1 



1-9* 5 (l 



i=0 

t-1 



qk 



A l : 



ir\lm— 1 



8 2m- 1 



l - qk 



i=0 v 



i /o2m— 1 



> + 



m— 1 



7=1 



k=l 



Taking into account that qt are roots of the Euler-Frobenius polynomial E 2m -2(q) and using 
theorem 2.1, for ^(xp) we have 

,, , , .1 r , Ml — 1 2(1) — ] -i W + i 

(<z* - 1) 5 



^2(^/3) 



/i 2m 

(2m- 1)! 



^ (2m-l)!E 2 ^ ^ - gfc + (-i)i^ 



t-i 



fc=l i=0 



m— 1 



7=1 V fc=i 



(4.10) 



Finally, using binomial formula and (2.6) from (4.10) we get 



, , , (h{3) 2m h 2m fi 2m - x 



m— 1 



9^ 1 """" n l - n N - t+l 



fc=i 



+ 



+ 



h 



2m 



2m-2 



{2m-l)\B j+l f3 2m - l " J h 2m 2 ^ (2m - ly.p 2171 ' 1 ^ 



X 



X 



"i-l J _ , / -, \i n N+i t-1 / m-1 

E * E =t^# A, ° : ' - E 1 + E * {-I + <T 7 

,fc=l i=l Wfe ' 7=1 V k=l 



(4.11) 



(2m -1)! ^ (j + l)!(2m- 1 -j)! (2m - 1)! j\(2m-l-j)\ 

"'»- J r i\i c N + , 

(q k - i) r 

Now we consider the equality (4.7). Using binomial formula, keeping in mind (1.5) and (3.3) 
when a — 0, for i])}_(xp) we have 



t-i 



t-i 



7=0 v 7 7=0 



(/l/3 — 

(2m- 1)! 



2m- 1 



(2m -1)! ^ j!(2m-l- J )! ^ ^ ( ^ + 



h 2mp2m-l 

' (2m- 1)! 



m— 1 



2t - 1 ^ . ,< - ,~ -« 



fe=i 



gfc 



(4.12) 
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Further using binomial formula and the equality (1.9) from (4.9) we obtain 

m-1 i i\i Jm-l-j m-1 / 1 \m+j m-l-j N 



E 



o 2(j + l)!(2m-l-j)! ^2(m + i)!(m-l-i)!^ (| 
Similarly, using binomial formula in the equality (4.5) for / m (x^) we have 



^ ^ C-yX 



(4.13) 



2m m-1 ( _iy+l x 2m - 1 -j 



+ E 



m— 1 



1-1-7')! ^ 



— l) m+j x m+1+j 



JV- j^2(m + l+j)\(m-l-j)l 



(4.14) 



(2m)! ^ 2 (j + l)!(2m 
Substituting equalities (4.11), (4.12), (4.13), (4.14) to (4.3) and after calculations we obtain 



h 2m 



2m-2 



yy (2m- l)\B j+1 fi 2m ~ l -J h 2m 2 ^ (2m - l)\p 2m ~ l ~i 
(2m- 1)! (j + l)!(2m- 1 - j)\ + (2m- 1)! j\(2m-l-j)\ * 

e (ck-y - y) + e * E - E (-^' (<s + * 



x 



7=0 



k=l 



+E 



3 1 \i„iV+i 

AT 



AT-7 



.m+j 

m + l+j ^ 



7=1 

•■Pm-l(^)- 



+ 



(4.15) 



^ 2(m + j)\(m-l-j)\ 
Keeping in mind the equality (4.13) and designations (4.4), (4.5), from (4.15) one can see that 



the difference g(xp) — f m (xp) is polynomial of degree 2m — 2 with respect to xp, i.e. 



2m-2 



g(xp) - fra(xp) = ^2 a i X p' t < f3 < N — t. 
j=0 



(4.16) 



Here 



bj as m < j < 2m — 2, 



bi + 



J 3 1 2j\(2m-j-l)\ 

k 2 ™ 



N 



(2m-l)! 



2m-j 
t-1 



7=0 



as 1 < j < m — 1, 



1 )2m-l £ (^-l^m-1 _ + 



7=0 



(4.17) 



m— 1 /2m— 1 , ,j jv+i t— 1 / \ 

E4( E '"1-^ ^ 2m ~ 1 + E 7 2 - 1 (g fc 7 + q N k -") 



k=l 



i=l 



where 



2(2m 



h 2m-j 

J ~ 3 \(2m- 3 -\)\ 



hiy. (2^ ~ T, o C 7 x ^ m ^ as 3 ~ 0> 



^2m- j + ( _ l) 2m-,-l |- ^-l^m-,-1 _ ^m-,-1) + 



2m — j 



7=0 
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m-1 /2m-j-l , w 7v+j * _1 \ 

+ E * ( E "V-i)^ fc AtQ2m ' 1 " + ( - 1)2m "^ 72m_H + 9 "~ 7 )) • 

On the other hand from (4.15) we obtain 



g{xp) ~ fm(xfs) = -P m -i(xp). (4.18) 
This equality takes placed if bj = 0, as m < j < 2m — 2 or 

•m—\ ( j / -\\% N+i t— 1 ^ 

E * {E =7^^ A ' y - i-'i'E' (a 7 + * } = 

t-i 

-D.7 + 1 



^-(-lyE^'^-y). J = 1,2,..., m-1. (4.19) 



^ + - 



From equalities (4.16) and (4.18) we will find unknown polynomial P m ^i(xp) of the system 
(1.8) - (1.9) 

m— 1 

p m -i(^) = -E a ^- ( 4 - 2 °) 

Later, applying to the sum A = l j [ol + 'h" 1 } formulas (2.4), (2.5), (2.6) and (3.2), we 

7=1 ^ ' 

obtain 

Substituting the equality (4.21) to the (4.19), after some simplifications we have 

m-1 f 3 P _ n N-t+l , t+i ) 

E * (-D* 1 E c? ( " E - ,„ _ ,U' a A '° p = 

k=l I p=0 i=0 Wfc > ) 

t-1 

rfj+i 



\ ~ (-1) J 'E ( c i h ~ X rt, ~ ^) > J = 1,2,..., m-1. (4.22) 

"7 7=0 

Multiplying the first equation of the system (4.22) by (— l) 2 C2t, adding with the second one, 
then, multiplying the first equation by {—VfClt 1 , the second by (— l) 2 C|t of the system (4.22), 
adding with the third and so on, continuing by this way, and also taking into account lemmas 
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3.2, 3.3, 3.4, 3.5 and binomial formula for unknown parameters dk we get following linear system 
of equations 



N 3 „t+l | ( -\\i„N-t+i .-■j-i n ; I 

E*E 



fc=l i=l ""' J ' /3=0 



j = l,2,...,m - 1. 

The last system is the system (4.2) for d k . This completes the proof of the theorem 4.1. 



4.2 Norm of the error functional of the optimal quadrature formula 

In this section square of norm of the error functional (1.2) of optimal quadrature formulas of 
the form (1.1) with the nodes (1.5) is calculated. 
Following is valid 

Theorem 4.2. Square of norm of the error functional (1.2) of optimal quadrature formulas 
of the form (1.2) with the nodes (1.5) on the space L^^O, 1) have following form 



*)|4 m) *(0,l) 



■1) 



m+l 



t-1 



1,2m p 0!,2m+l 



m— 1 2m 
fc=l i=0 



i„t+i JV-t+1 



9* 



i±2m 



where B a are Bernoulli numbers, Cp, (3 = 0,t — l are determined from the system (3.3), 



are the roots of the Euler-Frobenius polynomial E2 m -2(q), \qk\ < 1, Vp, P = 0,t — 1 are defined 
from (1.5), AH 2m is finite- difference of order i of t 2m , d& are determined from (4-2). 
Proof. Taking into account (4.3), we reduce the expression (1.7) to the form 



*)l4 m) *(o,i) =(-i) 



\m+l 



A? 



^2 Cp (fm(xp) + P m -i(xp)) ~ 



.(5=0 



[2m + l)\ 



(4.23) 



Applying the binomial formula to the equality (4.5), we obtain 

(1 -xp)-~ \ - I--' / 



fm(xp) 



2n, : ! ( — l) i+1 x^ m ~'~ :l 

(2m)! ^ 2(i + l)!(2m-i- 1)!' 



-E 



(4.24) 
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Putting (4.24) to the equality (4.23) and using equalities (4.17), (1.9), consecutively we have 



A? 



ii€(*)ii 2 = (-ir +i 

h 2m-j 



.13=0 



2m 



(2m)! (2m +1)! ^ (j + l)!(2m - j)\ 



m—l 



h 2m ~ j B. 



2m-j 



m—l 



t-1 



Z.( J + l)!(2m- J -l)! \ ( 1} ^ k ^ 7 j + 



3=0 

m—l /2m— j — 1 



t-1 



fe=i 



(_!)-+! [Z x + Z 2 ] , 



(4.25) 



N 



where Z\ = 2J C/3 ( 2 m)! — an< ^ ^ 2 * s ^ ne rema i n ing part in square brackets of the equality 

0=o ( m> ' 

(4.25). 

Keeping in mind symmetry of the nodes (1.5), and doing similarly calculations as in the 
proof of the lemma 3.3 when a = 2m, for Z\ we have 



Zi 



+ 



h 2m B 2m 2 ^ 2 B^h 2 ™^ 



+ E 



(2m + 1)! (2m)! (j + l)!(2m - j)\ 



+ 



2m-2 



+ £ 

3=0 



h 2m-j 



(j + l)!(2m-j-l)! 



t-i 



(_ l) 2m- J -l {Cph- 1 ^-*- 1 - P 2 ™-?- 1 ) + 



0=0 



m-1 (2m-j-l N+i , *v 

+E<M E fc ( } 



fe=i 



i=i 



(i - <? fc ) m 



t-i 

A * 2m-,-l + ( _ l) 2m- j ^ ^m-,-1 ^ + 



2/l 2m+ 

" (2m)! 



'm—l ( 2rrt AT+i 



fc=i I i=i 



c: + (-i) f g fc 

(i - 



0=0 

t-1 



+ 



0=0 



t-i 



+ E (^/T 1 ^ - P 2m ) ■ (4-26) 

/J=0 

Putting the equality (4.26) into (4.25) and taking into account the equality (4.19), we have 



\\i(x)\\ 2 = (-1) 



m+l 



h 2m B 



2m 



2h 2m+1 



(2m)! (2m)! 



m-l ( 2m N+i , y 

j2 d k \ e ~ % - + \.} 9k A i o 2m 



,k=l \i=l 



(i - q k y +i 
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E P 2m {4 + ) + E (^-^r - ? m ) 

13=1 ) 13=0 



(4.27) 



Hence, using formulas (2.4), (2.5) and lemma 3.2 when s = 0, after simplifications we will get 
the assertion of the theorem. 
Theorem 4.2 is proved. 

Remark. Should be noted, that when the nodes (1.5) are equal spaced, i.e. when in (1.5) 
Vo = 0, r/i = 1, ...,r] t -i — t — 1, from the theorems 4.1 and 4.2 we get the theorem 2.1 of [11] 
and the results of [31]. 



5 Numerical results 

For simplicity in this section we investigate the optimal quadrature formulas of the form (1.1) 
with the nodes 

x = i] h, x N = l-i] h, < i] < 1, Xf3 = h(3, (3 = 1, 2, iV - 1. (5.1) 
This means we consider optimal quadrature formulas of the form 

* AT-l 

J <p(x)dx = C <p(rj h) + E C P ip{hfi) + C N ip(\ - rjoh) (5.2) 



o 



(3=1 

with the error functional 



£(x) = e m (x) - (^C 5(x - r) h) + E C p 8(x - hp) + C N 5(x - (1 - r) h)) j . (5.3) 

Then for this case from theorems 4.1 and 4.2 when t — 1 we get following 

Corollary 5.1 The coefficients of optimal quadrature formulas of the form (5.2) with the 
error functional (5.3) in the Sobolev space L^^O, 1) are expressed by formulas 

(TO-l N \ 

, P = o, n, 

" > / / x\ (5 - 4) 

Ml+E d k (q% + q£-P) ) , [3 = 1,2,..., N-l, 
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where d k (k — l,m — 1) satisfy following system of m — 1 linear equations 



m—l 



"9fc + 



i)^ 1+ ; AI0J + fa-#)(i-^ y 



fe=i 



,i=l 



(<z* - i) m 



j = 1,2,.. .,771- l, 



j + l 2 ' - -' ( 5 - 5 ) 

/iere « r e the roots of the Euler-Frobenius polynomial E 2m -2(q) of degree 2m — 2, \q k \ < 1, Bj +1 
are Bernoulli numbers, < r/o < 1, A*7 J zs finite difference of order i of , A l (P = A l 7 J | 7= o. 

Corollary 5.2. Square of norm of the error functional (5.3) of optimal quadrature formulas 
of the form (5.2) on the space L 2 m ^(0, 1) have following form 



*)i4 m) *(o,i) 



■i) 



m+l 



2m 



2/l 2m+l 

' (2m)! 



m—l 



fc=i 



2m 

E 



"9* 



(2m)! 
A 4 2m + 



i=i 



(5.6) 



;i - <? fe ) l+1 g fc - i 

where B 2m is Bernoulli number, q k are the roots of the Euler-Frobenius polynomial E 2 m-2{q), 
\qk\ < 1, < 770 < I, AV m « /»ra*e difference of order i of-f 2m , A^O 2 " 1 = AV m | 7=0 . 

The aim of this section is by using the equalities (5.4), (5.5) and by choosing < r/ < 1 to 
investigate positiveness of optimal coefficients. 

As mentioned above, I.J.Schoenberg and S.Silliman in [25] showed, that in the case of equal 
spaced nodes among optimal coefficients it is appears negative coefficient starting from m = 7. 
From (5.4) and (5.5) in the case r/ = and iV — > oo taking instead of optimal coefficients Cp 
the coefficients NCp we get the results of [25]. 

By choosing the value of r/ = 0, 205 and using (5.4) and (5.5), we obtain optimal quadrature 
formulas of the form (5.2) with positive coefficients in the cases m = 2, 3, 14. 

Below, in the first subsection, we give the lists of the optimal coefficients in tables. In the 
second subsection, we compared some of the results of this paper with well known formulas. 
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5.1 The coefficients of optimal quadrature formulas of the form (5.2) in the case 

Vo = 0, 205 

Since \q k \ < 1 from (5.4) it is easy to see, that the optimal coefficients which located in the 
middle of the list of optimal coefficients are close to h. Therefore, and also taking into account 
symmetry of the optimal coefficients, we give the tables of optimal coefficients which satisfy 
the condition \Cp — h\ > 10~ 8 and are located in the left boundary layer in the list of optimal 
coefficients when N = 300, i]q = 0, 205 for m = 2,3, 14, 15. 



Table 1. The optimal coefficients for m = 2 



C = 0.00177612633884071008 
Ci = 0.00319454403039647954 
C 2 = 0.00337052181497334175 


C 3 = 0.00332336870971015344 
C 4 = 0.00333600334618604447 
C 5 = 0.00333261790554566866 


C 6 = 0.00333352503163128086 
C 7 = 0.00333328196792920787 
C 8 = 0.00333334709665188764 




Table 2. The optimal coefficients for m = 3 




C = 0.00175247906611553142 
Ci = 0.00324264900163471085 
C 2 = 0.00333892060971646866 


C 3 = 0.00333236955552006215 
C 4 = 0.00333368616893237390 
C 5 = 0.00333318408917128015 


C 6 = 0.00333339747877195358 
C 7 = 0.00333330571886293306 
C 8 = 0.00333334522322914572 




Table 3. The optimal coefficients for m = 4 




Co = 0.00174455012038852683 
Ci = 0.00326634746379700186 
C 2 = 0.00330774531081996009 
C 3 = 0.00335699483364953318 
C 4 = 0.00331942429918189420 
C 5 = 0.00334093116706873427 


C 6 = 0.00332924765792690917 
C 7 = 0.00333552260766573748 
C 8 = 0.00333216117670071227 
C 9 = 0.00333396080026924074 
Cio = 0.00333299745834218190 
Cn = 0.00333351312116039613 


C12 = 0.00333323709636443753 
C i3 = 0.00333338484710726915 
C14 = 0.00333330575901727716 
C15 = 0.00333334809332522744 


Table 4. The optimal coefficients for m = 5 


C = 0.00173901538774817543 
Ci = 0.00328716890432268279 
C 2 = 0.00326911260567167412 
C 3 = 0.00340205381465105221 
C 4 = 0.00328232898469587043 
C 5 = 0.00336635018632553123 
C 6 = 0.00331284808226362786 
C 7 = 0.00334587152301716133 


C 8 = 0.00332569334718594541 
Cg = 0.00333798181479271212 
Cio = 0.00333050638482074723 
Cn = 0.00333505224863080752 
C12 = 0.00333228820948218315 
Ci 3 = 0.00333396877152254017 
C14 = 0.00333294698744021598 
C15 = 0.00333356823085626853 


Ci 6 = 0.00333319051620654561 
C i7 = 0.00333342016578285375 
Cis = 0.00333328053942885252 
C19 = 0.00333336543188980042 
C 20 = 0.00333331381749472885 
C 2 i = 0.00333334519891226771 



Table 5. The optimal coefficients for m = 6 



Co 


= 0.00173498886058381026 


c 7 = 


= 0.00337747643535115270 


Cis 


= 0.00333497426498772165 


C 2 i 


= 0.00333347053324047546 


Ci 


= 0.00330520593286011352 


c 8 = 


= 0.00330385537900635393 


C14 


= 0.00333085188184676641 


c 22 


= 0.00333324260768654442 


c 2 


= 0.00322620390097277555 


Cg -- 


= 0.00335290464626403816 


Cl6 


= 0.00333224823220328272 


c 23 


= 0.00333339332712609729 


c 3 


= 0.00346783504823977116 


Cio 


= 0.00332037008745176751 


Cl 7 


= 0.00333405087626146100 


c 24 


= 0.00333329366147357165 


Ci 


= 0.00321158077766250064 


Cn 


= 0.00334191131490927107 


Cis 


= 0.00333285884589705156 


C 2 5 


= 0.00333335956698815529 


c 5 


= 0.00342625617820782061 


C i2 


= 0.003327659419014779572 


Cl 9 


= 0.00333364709595060379 


c 26 


= 0.00333331598590751689 


c 6 


= 0.00326815998330420462 


Ci 3 


— 0.00333708573208910388 


c 20 


= 0.00333312585271306499 


c 27 


= 0.00333334480459741380 
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Table 6. The optimal coefficients for m = 7 



Co 




0.00173205877707359127 


c 9 -- 


= 0.00339643943424144824 


Ol8 




0.00333069296984364163 


O27 


= 0.00333344250631211799 


Ci 




0.00332034204314634415 


C±o 


= 0.00328876219893682112 


O19 




0.00333518660336370406 


O28 


= 0.00333325670544635996 


c 2 




0.00318230520695443787 


On 


= 0.00336471020635444449 


C20 




0.00333203252906371234 


O29 


= 0.00333338711800696554 


c 3 




0.00355132503340265203 


Ol2 


= 0.00331127918647367579 


O21 




0.00333424636193460625 


O30 


= 0.00333329558217999618 


d 




0.00309977448873935734 


Ol3 


= 0.00334882331050388661 


C22 




0.00333269248328614613 


O31 


= 0.00333335983065088640 


c 5 




0.00354150500848728542 


Ol4 


= 0.00332245757624453130 


C 23 




0.00333378314247101608 


O32 


= 0.00333331473501844551 


c e 




0.00316818905922230260 


O15 


= 0.00334096810776062514 


C24 




0.00333301761482745647 


O33 


= 0.00333334638738364948 


c 7 




0.00345627307134070750 


Ol6 


= 0.00332797414827915069 


O25 




0.00333355493435702833 






Cs 




0.00324460285490503001 


O17 


= 0.00333709504135639172 


O26 




0.00333317779284220016 







Table 7. The optimal coefficients for m = 8 



Co 




0.00172995290066120350 


On 




0.00342842342134998057 


O22 




0.00333014109574784269 


O33 




0.00333343948477634317 


Oi 




0.00333261654096956311 


O12 




0.00326323297211348103 


O23 




0.00333567603784529139 


O34 


= 


0.00333325543170061418 


c 2 




0.00314030615358630034 


O13 




0.00338490036159433393 


O24 




0.00333161408331454461 


O35 




0.00333339050320497157 


c 3 




0.00364680181687277409 


C14 




0.00329544249146043423 


O25 




0.00333459504509005457 


O36 


= 


0.00333329137793257419 


C A 




0.00294599819613147818 


O15 




0.00336115859898192848 


O26 




0.00333240739717116319 


O37 


= 


0.00333336412325119880 


c 5 




0.00373067122924529534 


Ol6 




0.00331290610817136431 


O27 




0.00333401285268273812 


O38 


= 


0.00333331073745712191 


C 6 




0.00297727240591131934 


O17 




0.00334832702306404860 


O28 




0.00333283465264497077 


O39 


= 


0.00333334991582712020 


c 7 




0.00362534555471145307 


Ois 




0.00332232883088807872 


O29 




0.00333369930144181951 


O40 


= 


0.00333332116389597219 


Cs 




0.00310592534617643848 


Oig 




0.00334140963871753226 


O30 




0.00333306475936511862 








C 9 




0.00350552648624077583 


C20 




0.00332740619898615369 


O31 




0.00333353043240642612 








O10 = 


= 0.00320487233533179992 


o 2 i 




0.00333768315455876029 


O32 




0.00333318868772800817 










Table 8. The optimal coefficients for m 


= 9 












Co 




0.00172847784910360447 


012 




0.00314035977572450425 






0.00332610031017101433 


O36 


= 


0.00333306608365147515 


Oi 




0.00334216427951079832 


013 




0.00348084679101216452 


C25 




0.00333882817755420089 


O37 


= 


0.00333353635843483279 


c 2 




0.00310265417511838644 


Cn 




0.00322087398274449081 


C26 




0.00332915897604222655 


O38 


= 


0.00333317909856898933 


c 3 




0.00374645947596537845 


O15 




0.00341893779189990479 


C27 




0.00333650452957059022 


O39 




0.00333345050289666260 


c 4 




0.00275792168171570242 


016 




0.00326822729185132696 


C28 




0.00333092422572131367 


O40 




0.00333324432158157594 


c 5 




0.00400138156396844176 


017 




0.00338282511219261974 


C29 




0.00333516349327706089 


O41 




0.00333340095406794224 


Co 




0.00266165221372250548 


Ois 




0.00329572154562808534 


C30 




0.00333194299096395421 


O42 




0.00333328196299562843 


c 7 




0.00394079887384133675 


O19 




0.00336191228831355167 


C31 




0.00333438955334958306 


O43 




0.00333337235851711259 


Cs 




0.00282156182791894007 


O20 




0.00331161983320458289 


C32 




0.00333253094060154966 


O44 




0.00333330368655589561 


Co 




0.00374603753952685103 


O21 




0.00334982981298158293 


C33 




0.00333394289767582404 


O45 




0.00333335585549280241 


O10 = 


= 0.00300894501757551301 


O22 




0.00332080076136125273 


C34 




0.00333287025750894230 


O46 




0.00333331622362661400 


On = 


= 0.00358457992357594231 


O23 




0.00334285432190737266 


^35 




0.00333368512427696033 


O47 




0.00333334633129049714 




Table 9. The optimal coefficients for m 


= 10 










Co 




0.00172749278862645234 


Cn 




0.00304088010987428129 


028 




0.00332404459372417715 


O42 




0.00333304181986959795 


Oi 




0.00334915954552527705 


Ois 




0.00356281566915249473 


029 




0.00334058737008624209 


O43 




0.00333356098973720558 


c 2 




0.00307141154907661477 


C16 




0.00315360445472180508 


O30 




0.00332766830567799599 


O44 




0.00333315554587373269 


c 3 




0.00384087716305041600 


017 




0.00347393491400228088 


031 




0.00333775742286366819 


O45 




0.00333347217581801177 


Ci 




0.00255327063613243166 


Ois 




0.00322341661691745051 


032 




0.00332987835336878476 


O46 




0.00333322490478845213 


o 5 




0.00434013033335869961 


O19 




0.00341922625493328018 


033 




0.00333603148905263082 


O47 




0.00333341801022001505 


Co 




0.00221013903097713541 


O20 




0.00326623022080639272 


034 




0.00333122621755866373 


O48 




0.00333326720521970105 


c 7 




0.00444946400666106578 


O21 




0.00338574919280914013 


035 




0.00333497887821597101 


O49 




0.00333338497584858724 


Cs 




0.00231492992409806705 


O22 




0.00329239374555605829 


036 




0.00333204825076180440 


O50 




0.00333329300329463880 


Co 




0.00420985817056583358 


O23 




0.00336530759722583348 


037 




0.00333433691406739092 


O51 




0.00333336482893367433 


O10 = 


= 0.00260680747387077026 


O24 




0.00330836190455046297 


038 




0.00333254959049473769 


O52 




0.00333330873695621772 


On = 


= 0.00392158564879464832 


O25 




0.00335283526564426801 


039 




0.00333394539453196409 


O53 




0.00333335254178659839 


012 - 


= 0.00286380539858292753 


O26 




0.00331810309402904888 


O40 




0.00333285534630910409 


O54 




0.00333331833256004469 


013 - 


= 0.00370486410231191268 


O27 




0.00334522746373694056 


041 




0.00333370661560394268 


O55 




0.00333334504813390320 



23 



Table 10. The optimal coefficients for m = 11 



Co = 


0.00172689334490066565 


C\q 







00290248747287229047 


C32 







00332138043174862451 


Ci8 







00333300511542151927 


Ci = 


0.00335378561374918185 


C±7 







00367878243460793137 


C33 







00334288089094969269 


Cio 


= 





00333359550162011932 


C 2 = 


0.00304832508611045170 


Cis 







00305674037361827340 


C3i 







00332570708681864833 


C50 







00333312392303583896 


c 3 = 


0.00391935697738649923 


C±o 







00355459988011267817 


C35 







00333942490198512218 


C51 


= 





00333350060250695466 


Ci = 


0.00236091789799528996 


C20 







00315642482727045543 


C36 







00332846761136009168 


^52 


= 





00333319972492547481 


c 5 = 


0.00470103221229688833 


C21 







00347472666712902671 


C37 







00333721989298164747 


C53 


= 





00333344005477077078 


c 6 = 


0.00166606164689116412 


C22 







00322035045788260487 


C38 







00333022889296381827 


C^i 


= 





00333324808820998967 


c 7 = 


0.00513715438246413703 


C23 







00342360164155868860 


C30 







00333581304546978635 


C55 


= 





00333340142397405117 


C 8 = 


0.00155687727226540244 


C24 







00326121936040494949 


Cio 







00333135263111543830 


C56 


= 





00333327894505666322 


c 9 = 


0.00496522549139354527 


C25 







00339094085206410810 


Cn 







00333491544485158395 


C57 


= 





00333337677667015449 


Cio = 


= 0.00190545757664787945 


C20 







00328731578129565660 


Ci 2 







00333206960132156587 


C58 


= 





00333329863240825021 


Cn = 


= 0.00454262920324724638 


C27 







00337009186182822834 


C43 







00333434275558700776 


Csg 


= 





00333336105114438844 


C12 - 


= 0.00233068337371350425 


C28 







00330397127017939317 


C44 







00333252704428838619 


Ceo 


= 





00333331119337459038 


C\3 - 


= 0.00415339159536806206 


C20 







00335678702055923429 


C45 







00333397736710399180 


Coi 


= 





00333335101791066748 


C14 - 


= 0.00266841337297767427 


C30 







00331459921621225224 


Cio 







00333281890305696659 


C&2 


= 





00333331920754931007 


Cl5 ~ 


= 0.00386950335147916675 


C31 







00334829754279374992 


C47 







00333374424114699843 


C&3 


= 





00333334461648534207 


Table 11. The optimal coefficients for m = 


12 












Co = 


0.00172660093000204434 


Cis 







00275760059774337822 


C 36 







00331906522131345626 


C S 4 


= 





00333298282252256406 


Ci = 


0.00335622016635424541 


C±o 







00380333655344864095 


C37 







00334494624845792310 


C55 


= 





00333361861605357889 


C 2 = 


0.00303488707201474333 


C20 







00295003782640532812 


C38 







00332388150649666917 


C56 







00333310114012289132 


c 3 = 


0.00397016939076339731 


C21 







00364570570477646406 


C3Q 







00334102623365402826 


C57 


= 





00333352231669313408 


Ci = 


0.00222172391357403828 


C22 







00307887414398914889 


Cio 







00332707203566076212 


C$8 







00333317951871399489 


c 5 = 


0.00499378828442308939 


C23 







00354055563717261057 


Cn 







00333842944024015643 


C50 


= 





00333345852391701871 


C 6 = 


0.00117134786383522544 


C24 







00316461152396595462 


Ci 2 







00332918558244974745 


Coo 


= 





00333323144001324420 


C 7 = 


0.00583516570203082313 


C25 







00347069038125884784 


Ci3 







00333670921153828293 


Co\ 


= 





00333341626487956293 


c& = 


0.00070537581513706783 


C20 







00322151976711594156 


Ca 







00333058568684464716 


C02 


= 





00333326583488148334 


C 9 = 


0.00589276164637635086 


C27 







00342434860282943662 


C45 







00333556965838687970 


Cq3 


= 





00333338827070287737 


Cio = 


= 0.00097691145372704574 


C 2S 







00325925036543305965 


Cio 







00333151317572930858 


Coi 


= 





00333328861949931687 


Cn = 


= 0.00541764827809579119 


C20 







00339363262285636622 


C47 







00333481476987879354 


Con 


= 





00333336972617429080 


Cl2 - 


= 0.00154000801804920969 


C30 







00328425393283165749 


Ci8 







00333212758377888206 


Ceo 


= 





00333330371299542040 


Cl3 = 


= 0.00484769113306215419 


C31 







00337328012014849898 


Cio 







0033343 1469968656027 


Ce7 


= 





00333335744149293259 


C14 - 


= 0.00207043279270397406 


C32 







00330081999773997306 


C50 







00333253459373417958 


Cos 


= 





00333331371156666169 


Cl5 - 


= 0.00437783417724906036 


C33 







00335979634546885836 


C51 







00333398343198778386 


Ceo 


= 





00333334930359957163 


C\q - 


= 0.00247418267818456252 


C3i 







00331179481824981668 


C52 







00333280421438103343 


C70 


= 





00333332033504383217 


C\j - 


= 0.00403747731479976406 


C35 







00335086371649026079 


C53 







00333376398622498520 


C71 


= 





00333334391271425501 


Table 12. The optimal 


coefficients for m = 


13 












Co = 


0.00172655548945567044 


C20 




0.00277237192137831871 


Cio 




0.00332071149652304087 


Ceo 




0.00333305114925564172 


Ci = 


0.00335662901376447258 


C21 




0.00379816301931711370 


Cn 




0.00334377075895046149 


Cei 




0.00333356668071934869 


C 2 = 


0.00303238425807657517 


C22 




0.00294841399870425801 


Ci 2 




0.00332470227406376349 


Ce 2 




0.00333314037062970202 


c 3 = 


0.00398073517633032120 


C23 




.00365 193776947454 155 


Ci3 




0.00334047064592444503 


Ce3 




0.00333349290060179606 


Ci = 


0.00218920482037733280 


C 2 i 




0.00306969888569782702 


Ca 




0.00332743125080413949 


Cei 




0.00333320138185937509 


c 5 = 


0.00507104524959419857 


C25 




0.00355143697213385530 


Ci5 




0.00333821396235558649 


Cos 




0.00333344244838961507 


Co = 


0.00102336156051060256 


C20 







00315292244059992812 


Cio 




0.00332929737855920021 


Cee 







00333324310246929006 


C 7 = 


0.00607199798241084696 


C27 




0.00348255132276115587 


Ci7 




0.00333667079868275828 


Ce7 




0.00333340794823161848 


c 8 = 


0.00037867881870999653 


C28 




0.00320992302093593156 


Ci8 




0.00333057347222005968 


C 6 8 




0.00333327163179578226 


Co = 


0.00629239692582258769 


C20 




0.00343539498005851424 


Cio 




0.00333561555424564661 


Ceo 




0.00333338435638400590 


Cio - 


= 0.00053211065682692867 


C30 




0.00324892979162245131 


C50 




0.00333144608900867865 


C70 




0.00333329114067812168 


Cn = 


= 0.00587832107921238828 


C31 




0.00340313244846696006 


C51 




0.00333489395837593862 


C71 




0.00333336822384125740 


Cl2 - 


= 0.00108774569344314697 


C32 




0.00327561242680296621 


C52 




0.00333204280063615644 


C72 




0.00333330448121352504 



24 



Table 12. Continuation 



Cl3 


= 0.00527471154107237576 


C33 


= 0.00338106562394459128 


C53 


= 0.00333440051769216630 


C73 


= 0.00333335719210945782 


C14 


= 0.00167844383373434332 


C34 


= 0.00329386137938002601 


C54 


= 0.00333245084313764205 


C74 


= 


00333331360371894311 


C15 


= 0.00473043759244898944 


C35 


= 0.00336597430528210186 


C55 


= 0.00333406309378899346 


C75 


= 0.00333334964840674512 


Cl6 


= 0.00216162520663729902 


C36 


= 0.00330634125957574669 


C56 


= 0.00333272987017358315 


C70 


= 0.00333331984185726252 


C17 


= 0.00431159258197410437 


C37 


= 0.00335565407117066211 


C57 


= 0.00333383235706368672 


C77 


= 0.00333334448988298376 


Cis 


= 0.00251908554798301532 


C38 


= 0.00331487552014629517 


C58 


= 0.00333292067403627449 








C19 


= 0.00400965449172913550 


C39 


= 0.00334859674230775176 


C59 


= 0.00333367457501187563 








Table 13. The optimal coefficients for m = 


14 










Co -- 


= 0.00172671049778565797 


C23 


= 


00316025452220934915 


C4Q 


= 


00333638193499156475 


C09 


= 


00333328063425493357 


Ci = 


= 0.00335516395507917580 


C24 


= 


00347924238639198785 


C47 


= 


00333077781143822645 


C70 


= 


00333337750881443832 


c 2 -- 


= 0.00304193627584372803 


C25 


= 


00321053936117372331 


C48 


= 


00333547552503919473 


C71 


= 


00333329630283508029 


c 3 -- 


= 0.00393776327125170017 


C20 


= 


00343655097788938971 


C49 


= 


00333153762035049978 


C72 


= 


00333336437448481488 


Ci -- 


= 0.00233007712879415692 


C27 


= 


00324664320823834275 


Cso 


= 


00333483860682981228 


C73 


= 


00333330731280618138 


c 5 = 


= 0.00471527133757507641 


C28 


= 


00340610011591843955 


C51 


= 


00333207152369661671 


C74 


= 


00333335514527654166 


c 6 -- 


= 0.00174497854628782157 


C29 


= 


00327227842278880512 


C52 


= 


00333439105698931688 


C75 


= 


00333331504927383940 


c 7 -. 


= 0.00485675294324682823 


C30 


= 


00338454689233366739 


C53 


= 


00333244668671512832 


C70 


= 


00333334866011226255 


Cs -- 


= 0.00212747598467839435 


C31 


= 


00329038332124877627 


C54 


= 


00333407657300585164 


C77 


= 


00333332048552169609 


Cg -- 


= 0.00408453051443764890 


C32 


= 


00336934817236508896 


C55 


= 


00333271030587360587 


C78 


= 


00333334410312812365 


Cio 


= 0.00303849317411607795 


C33 


= 


00330313676378294222 


C56 


= 


00333385559179154898 


C79 


= 


00333332430545503779 


Cn 


= 0.00326150825980108299 


C34 


= 


00335864987188439309 


C57 


= 


00333289554541894104 


Cgo 


= 


00333334090103430680 


C12 


= 0.00365086987677036012 


C35 


= 


00331210917522999367 


C58 


= 


00333370031303950707 


C$i 


= 


00333332698963926159 


C13 


= 0.00288027997265296082 


C36 


= 


00335112604866268579 


C59 


= 


00333302570925529457 


C&2 


= 


00333333865099250313 


C14 


= 0.00383972630105596425 


C37 


= 


00331841762354822903 


Ceo 


= 


00333359120201437040 


C&3 


= 


00333332887575723863 


C15 


= 0.00282694633373975373 


C38 


= 


00334583703432045650 


C01 


= 


00333311717256726671 


C%4 


= 


00333333706993697896 


Cl6 


= 0.00380915983064439615 


C39 


= 


00332285171576454862 


C&2 


= 


00333351453205997142 


^85 


= 


00333333020109139106 


C17 


= 0.00290295426940039930 


C40 


= 


00334211980668418472 


Ce3 


= 


00333318144185430803 


Cso 


= 


00333333595896358890 


Cl8 


= 0.00371316819510743415 


C41 


= 


0033259678971 1823063 


C&4 


= 


00333346065774556324 


Cg7 


= 


00333333113237492521 


C19 


= 0.00300351110208944757 


C42 


= 


00333950753291683119 


Cea 


= 


00333322660248814852 


Cgg 


= 


00333333517830670844 


C20 


= 0.00361660703987917620 


C43 


= 


00332815772174759486 


Ce,6 


= 


00333342280143380223 








C21 


= 0.00309185019337987811 


C44 


= 


00333767185593130795 


C07 


= 


00333325833588546756 








C22 


= 0.00353813630227038258 


C45 


= 


00332969651611532286 


Cos 


= 


00333339620062519364 









Table 14. The optimal coefficients for m = 15 



Co = 


0.00172702945455083193 


C20 







00591821410669587006 


C52 







00336928787611088775 


C78 







00333383067161036060 


Ci = 


0.00335196263611974402 


C27 







00113739105014957971 


C53 







00330283683648456223 


C79 







00333291149314649690 


C 2 = 


0.00306452563504523993 


C 2S 







00519802682335201819 


C54 







00335920033583227050 


C 8 o 







00333369113636488924 


C 3 = 


0.00382735744574236094 


C29 







00175042759986189166 


C55 







00331139305001559083 


C 8 i 







00333302984634137357 


C4 = 


0.00272471104016856072 


C30 







00467672632259750788 


C56 







00335194298912504867 


C 8 2 







00333359074969419812 


c 5 = 


0.00362523740203732478 


C31 







00219339925074939015 


C57 







00331754870221566789 


C83 







00333311499388354665 


Co = 


0.00416828584678143709 


C32 







00430050875356561753 


C58 







00334672179022929618 


C 8 4 







00333351852791527129 


C7 = 


0.00038025393759564147 


C33 







00251280401810681106 


C59 







00332197730137501015 


C 8 5 







00333317625208569821 


Cg = 


0.00918605541537786417 


C34 







00402940935773975411 


Ceo 







00334296547103844503 


C 8 6 







00333346656897177222 


C9 = 


-0.0056482216521837456 


C35 







00274285983082156931 


Coi 







00332516339541746467 


C 8 7 







0033332203234380 1599 


C10 - 


= 0.01503979792439844429 


C30 







00383420982729031424 


C02 







00334026303961781421 


C 8 8 







00333342918784036036 


Cn = 


= -0.0102513216994999018 


C37 







00290846808379631022 


Ce3 







003327455586219771 1 1 


C89 







00333325202995283771 


C12 = 


0.01780870154874396918 


C 38 







00369371677971561177 


C&4 







00333831881318149643 


C90 







00333340229451091797 


C\3 = 


-0.0111516599198095529 


C39 







00302764878602438014 


Cos 







00332910467075450735 


C91 







00333327484075914716 


C14 = 


= 0.01717584313564491950 


C40 







00359261903731864099 


Coo 







00333692006675210665 


C92 







00333338294648316812 


C15 = 


= -0.0094596822947696162 


C41 







00311340503994861123 


C07 







00333029108147133945 


C93 







00333329125167301758 


ClG = 


0.01487372890505018792 


C42 







00351987755666928677 


Cos 







00333591375827837438 


C94 







00333336902681665769 


C\7 - 


-0.OO68972585187668098 


C43 







00317510611508734637 


C09 







00333114462799414914 


C95 







00333330305827423077 



25 



Table 14. Continuation 



Cl8 


= 0.01228998089227584958 


C44 


= 0.00346754171514042525 


C70 


= 0.00333518978378703653 


Cqq — 


0.00333335901250834411 


C19 


= -0.0044379368869592209 


C45 


= 0.00321949791378515819 


C71 


= 0.00333175870004468364 


C<)7 — 


0.00333331155236741126 


C 20 


= 0.01003287196894858974 


C46 


= 0.00342988831851899784 


C72 


= 0.00333466893051466329 


^98 — 


0.00333335180785510696 


C21 


= -0.0024157213542609094 


C47 


= 0.00325143564488660772 


C73 


= 0.00333220048555005048 


C99 = 


0.00333331766332322668 


C22 


= 0.00825049420457872496 


C48 


= 0.00340279870089829129 


C74 


= 0.00333429420994935355 


Cioo - 


= 0.0033333466245687503 


C23 


= -0.0008623772730946671 


C49 


= 0.00327441304627985970 


^75 


= 0.00333251832192158444 


Caoi - 


= 0.0033333220597643263 


C24 


= 0.00690738675157999725 


C50 


= 0.003383309304493717340 


C76 


= 0.00333402462250633817 






C25 


= 0.OO029252137237138985 


C51 


= 0.00329094390700659850 


C77 


= 0.00333274698483095824 







The tables 1-13 confirm positiveness of the optimal coefficients of optimal quadrature for- 
mulas of the form (5.2) with the nodes (5.1) when t]q = 0, 205. 

5.2 Comparison some of the results with well known formulas 

In this section we compare some of the results of this work with the example (e) of [12]. From 
the example (e) of [12] when [a, b] = [0, 1] and for the equidistant nodes 

xp = hp, /3 = 0,1, ...,7V, h = ^ 

of the interval [0, 1] we obtain following particular formulas which were already obtained by 
Sard in [14, 18]. 

If N — 2 the optimal quadrature formula of the type (1.1) is 
1 

J V (x)dx = ^ [3y?(0) + l(V(l/2) + 3y?(l)] + R[ip] (5.7) 


and for the remainder we have following estimate 

\R[(p]\ < |M L 2 2) (°> 1)11 • 0.01398. (5.8) 

For N = 3 the optimal quadrature formula is 
1 

J V {x)dx = -i [4^(0) + 1M1/3) + 11^(2/3) + 4y?(l)] + R[<p] (5.9) 


and for the remainder R[<f\ following estimate is valid 

\R[<p]\ < \\ip\Lf\0, 1)|| • 0.00586. (5.10) 
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When N = 4 the optimal quadrature formula is 

1 

J <p(x)dx = ^ [11(^(0) + 32^(1/4) + 26^(1/2) + 32y?(3/4) + 11^(1)] + R[(p] (5.11) 
o 

and for the remainder R[<p\ we get 

|%>]| < ||^|4 2) (0, 1)|| • 0.00305. (5.12) 

For iV = 5 we obtain following optimal formula 

i 

J V {x)dx = ^[15^(0) + 43^(1/5) + 37y?(2/5) + 
o 

+37y?(3/5) + 43y?(4/5) + 15^(1)] + R[<p] (5.13) 

with 

\R[<p]\ < |bl4 2) (0, 1)|| • 0.00188. (5.14) 
Now we consider the quadrature formulas of the form (5.2). 

From (5.2), (5.4), (5.5), (5.6) when m = 2 and r] = for N = 2,3,4,5 we get the optimal 
quadrature formulas (5. 7)- (5. 14), respectively. 

From (5.2), (5.4), (5.5), (5.6) for m = 2 and r] = 0.205 when iV = 2, 3, 4, 5 we get following 
optimal quadrature formulas. 

For N = 2 the optimal quadrature formula is 
i 

J ip{x)dx = C ip(r] h) + CMl/2) + C 2V 9(1 - rjoh) + R[ip), (5.15) 
o 

where 

C = 0.27075812274368231046, 

d = 0.45848375451263537906, 

C 2 = 0.27075812274368231046 
and for the remainder R[(p] following is valid 

\R[(p]\ = \(£,(p)\ < |kl4 2) (0, 1)|| 2 -0.00694814. (5.16) 



27 



For N = 3 the optimal quadrature formula is 
i 

J ip(x)dx = C o ip(r] h) + C^(l/3) + C 2 <p(2/3) + C 3 <p(l - r] h) + R[<p], (5.17) 
o 

where 

C = 0.17683465959328028293, 

Ci = 0.32316534040671971706, 

C 2 = 0.32316534040671971706, 

C 3 = 0.17683465959328028293 
and 

\R[<p]\ = \(£,<p)\ < ||v?|4 2) (0,l)|| 2 • 0.00340515. (5.18) 

For N = 4 the optimal formula is 
i 

J <p(x)dx = C <p(voh) + C^l/4) + C 2V (l/2) + C 3 y?(3/4) + C A <p(l - Voh) + R[<p], (5.19) 
o 

where 

C = 0.13336566440349175557, 

Ci = 0.23884578079534432589, 

C 2 = 0.25557710960232783705, 

C 3 = 0.23884578079534432589, 

C 4 = 0.13336566440349175557 
and 

\R[(p]\ = \(£,<p)\ < ||v?|4 2) (0, 1)| | 2 - 0.0020343. (5.20) 

For N = 5 the optimal quadrature formula is 
i 

J <p(x)dx = C oV (r] h) + C^(l/5) + C 2V ?(2/5) + C 3 y?(3/5) + 
o 

+C 4 y?(4/5) + C 5 <p(l - Voh) + R[<p], (5.21) 
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C = 0.10653409090909090909, 

d = 0.19183238636363636363, 

C 2 = 0.20163352272727272727, 

C 3 = 0.20163352272727272727, 

C 4 = 0.19183238636363636363, 

C 5 = 0.10653409090909090909 
and for the remainder following estimate is valid 

\R[(p]\ = \(£,<p)\ < ||v?|4 2) (0,l)|| 2 - 0.0013408. (5.22) 

Thus, hence clear that the errors (5.16), (5.18), (5.20), (5.22) of the quadrature formulas 
(5.15), (5.17), (5.19), (5.21) are less than the errors (5.8), (5.10), (5.12), (5.14) of the quadrature 
formulas (5.7), (5.9), (5.11), (5.13), respectively. 
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